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I. INTRODUCTION 



The bulk viscosity of matter in neutron stars has recently received considerable attention in 
connection with damping of neutron star pulsations and gravitational radiation driven instabilities, 
especially in the damping of r-modes ,1,]. If the r-modes are unstable, i.e. if the damping time- 
scales due to viscous processes in neutron star matter are longer than the gravitational radiation 
driving time-scale, a rapidly rotating neutron star could emit a significant fraction of its rotational 
energy and angular momentum asgravitational waves, which could be detectable. 

It has been shown ^ Q, 0, Q, IJthat the bulk viscosity in a neutron star is caused by energy 
dissipation associated with nonequilibrium weak interaction reactions in a pulsating dense matter. 
Strong interaction processes do not play a role, because the strong interaction equilibrium is reached 
so fast that these processes can be considered to be in thermal equilibrium compared to the typical 
pulsation time-scales of 10~^ — 10~^ s. The relevant nonequilibrium weak interaction reaction rates 
and the magnitude of the bulk viscosity depend on the density and the composition of neutron 
star matter. 

In a relatively low density neutron star composed mainly of neutrons n with admixtures of 
protons p, electrons e, and muons fi, the bulk viscosity is mainly determined by the reactions of 
the nonequilibrium modified Urea process, 

N + n^N + p + l+Vi, N + p+ l^N + n+ui, (1) 

with N = n,p and I = e,fi. The bulk viscosity was studied by Sawyer [31 for npe matter and by 
Haensel et al. Q| for npefi matter. For the process in Eq. the relaxation time, i.e. the time it 
takes to restore equilibrium in case of a perturbation, is strongly temperature dependent, namely 

At densities of a few times saturation density, no {no ~ 0.15 fm~'^), the direct Urea process, 

n ^ p + I + Vi, p + I ^ n + vi, (2) 

may also be allowed depending on whether or not the proton fraction exceeds the Urea limit of 
n 

Xp ~ 0.11 1^]. The contribution of this process to the bulk viscosity was computed by Haensel 
and Schaeffer for npe matter and by Haensel et al. Q] for npefi matter. Compared to the 
modified Urea process a smaller number of particles is involved which leads to a weaker temperature 
dependence, r^^ oc T^. As a consequence at typical neutron star temperatures, T ~ 10^ — 10^'^ K, 
its contribution to the bulk viscosity is typically 4 — 6 orders of magnitude larger than that from the 
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modified Urea proeess. The largest difference compared witli tliat from the modified Urea process 
is reached at the low temperatures. 

At about the same densities hyperons may appear in the neutron star core, first the S~ and A 
hyperons, followed by H*^, H^, and S+ at higher densities. Here we will restrict ourselves to the 

and A hyperons. In addition to the semileptonic hyperon processes, weak nonleptonic hyperon 
processes also occur, specifically the processes 

n + n p + S^, (3) 



p + n <^ p + A, (4) 

and 

n + n <-> n + A. (5) 

At low temperatures these processes contribute more efficiently to the bulk viscosity than the direct 
Urea process and the semileptonic hyperon ones, because they contain no neutrino phase space 
factor; at typical neutron star temperatures, T < 10^'' K, the phase space of neutrinos is almost 
negligible compared to that of baryons. Hence, for the weak nonleptonic hyperon processes of 
Eqs. the temperature dependence of the inverse relaxation time is r^^ oc T^. 

Historically the first semi-quantitative calculation of bulk viscosity in neutron matter was carried 
out by Jones 10]. In this and all later works the weak nonleptonic process was calculated using a 
baryon current-current interaction, i.e. a contact W exchange. 

More recently the contribution to the bulk viscosity from the various weak nonleptonic hyperon 
processes has been reconsidered by several authors using a modern equation of state (EoS). Haensel 
et al. iJl studied the bulk viscosity for the nnpT,^ process within the nonrelativistic limit, and they 
found the bulk viscosity to be several orders of magnitude larger than that of the direct and the 
modified Urea processes. Also applying the contact interaction of W exchange, Lindblom and 
Owen '(Ji calculated the contribution of the pnpA process in addition to that of the nnpT," process. 

One may question the validity of the W exchange process. First there is no W exchange contribu- 
tion to the nnnA process and therefore it has not been considered in the above approaches P,|6, 10 1. 



On the other hand it is well known that in nuclear physics experiments on weak A decay in large 
hypernuclei the nnnA process is found equally important as the pnpA process. Secondly, Jones 
noticed already several orders of magnitude difference between Tnnnk and the theoretical W ex- 
change based value of T^npT,- iSl]- This is interesting, because the rate for the nnnA process must 



depend on a very wide class of possible weak hadronic processes. Because most of these processes 
also contribute to nnpT,^ process, the r^npT,- reality is probably of the same order of magnitude 

However, in nuclear physics it is customary to model the weak nonleptonic NY NN process 
in terms of meson exchange (pion and kaon) with one phenomenological weak vrA^y or KNN 
vertex. In this approach one can describe both the observed weak decay rate of the A as well as 



branching ratios in hypernuclei quite well 



Therefore, in the following we will calculate the 



bulk viscosity using the meson exchange picture to describe the hyperon processes in Eqs. (jSHHl)- 

In practice the rates also depend on the details of the equation of state, e.g. the hyperon frac- 
tions as a function of the density. A variety of model equations of state have been constructed 
with widely varying properties during the last decades, some of these based upon a microscopic 
free space A^A^ and NY interactions, while others are phenomenological parametrizations of the 
energy density to higher densities. In particular the relativistic mean field approaches, as used in 



Ref. Ul, and Ref. £1, although microscopic in nature, do not start from a realistic nucleon-nucleon 

n 

interaction. Here we employ the effective EoS based on the work of Balberg and Gal 12]. To 
get some idea about the sensitivity of the bulk viscosity to the details of the EoS we consider two 
different parameter sets for the density dependence of the multi-body potential energy. 

Superfluidity can also influence the bulk viscosity. Below the critical temperature, superflu- 
idity suppresses the reaction rates by roughly a factor exp(— A/T) with gap energy A. It 
generally leads to a smaller bulk viscosity. Above the critical temperature, it has no influence on 
the reaction rates. Many properties of superfluid matter such as the gaps and the critical temper- 
atures are still known with large uncertainties, e.g. Tc ~ 10® — lO^'^-fC. Hence, the bulk viscosity is 
only considered here for nonsuperfluid matter. For recent papers about the effects of superfluidity 
we refer to Haensel et al. and Lindblom and Owen Q]. 

The main goal in this study is threefold (i) to compute the viscosity using the more realistic 
one-pion exchange (OPE) instead of the contact W exchange description in the nonleptonic weak 
hyperon processes, (ii) to include the weak neutral current nnnA process, and (iii) to examine the 
sensitivity of the bulk viscosity to the EoS. The starting point of our derivation of the bulk viscosity 
for the nnpTi~ , pnpA, and nnnA processes is the finite temperature Green's function formalism in 
the quasi-particle approximation (QPA). The coupling constants for these processes using OPE can 
be verified from pionic hyperon decay experiments. The application of these coupling constants 
gives the correct order of magnitude for the rate of nonmesonic hyperon decay in large hypernuclei. 
Covariant expressions for the matrix elements are derived. 
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FIG. 1: The particle fraction as a function of the baryon number density (Balberg and Gal case 1 with 
7 = <5 = 2). 



The plan of this paper is as follows. The EoS is discussed in Section ^ The bulk viscosity 
for neutron star matter is discussed in section IIIII Section IIVI is devoted to the collision rate of 
the weak nonleptonic hyperon process as in Eqs. H3I5|) . The results are presented and discussed in 
section 13 Finally the conclusion is given in section IVTl 



II. EQUATION OF STATE 



121. It is a ffen- 



In this paper we use the equation of state constructed by Balberg and Gal 
eralization of the (npe) EoS of Lattimer and Swesty _13|. The energy density is parametrized in 
terms of the densities of the constituents. The parameters for the nucleon-nucleon interactions are 
fitted to the binding energy, symmetry energy, and incompressibility of saturated symmetric nu- 
clear matter. The values of the parameters for hyperon-nucleon and hyperon-hyperon interactions 
mainly rely on experimental data from hypernuclei. The equilibrium fractions are obtained in the 
standard manner, by assuming a given nucleon density n and solving the imposed conditions of 
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FIG. 2: The particle fraction as a function of tlie baryon number density (Balberg and Gal case 2 with 
7 = <5 = 5/3). 



baryon conservation, charge neutrality and weak equilibrium 



UB = Tin + Up + ns- + riA, (6) 

Up = Ue + ns-, (7) 

fJ-p = fJ-n - IJ-e, (8) 

/^e = ^J'^l, (9) 

AiE- = Ain + fJ-e, (10) 

AiA = /^n- (11) 



To obtain some idea about the sensitivity of the resulting bulk viscosity to details of the EoS we 
use two parametrizations which differ in the density dependence of the contributions of the nucleons 
{5) and hyperons (7) to the energy density, specifically given by the parameter sets 7 = 5 = 2 
and 7 = 5 = 5/3, respectively. The difference between the two parametrizations of the equation 
of state has no influence on the particle fractions in the nuclear sector. For densities below the 



appearance density of the hyperons the neutrons are abundantly present (x„ > 0.9). The leptons 
are present because of charge neutrahty to compensate the positively charged protons, although 
they are expensive in terms of energy density. The appearance of the hyperons create the possibility 
to lower the neutron excess without lepton formation, whereas the negatively charged hyperons 
allow charge neutrality to be maintained within the baryon sector. Therefore, the appearance of 
the hyperons is accompanied by a strong deleptonization. 

The first hyperon to appear is the S", followed by the lower-mass A. This can be explained 
by the higher chemical potential of the T,~ due to its negative charge, which compensates the 
mass difference of about 80 MeV. However, the growth of the T,~ number density fraction is soon 
hindered by charge dependent forces, which disfavor an excess of S~'s over S+'s and a joint excess 
of S~'s and neutrons. The A is not affected by charge-dependent forces and its fraction continues 
to increase. For the EoS with the smaller values for 5 = 7 = 5/3 the onset of the various hyperons 
occur at higher densities and the corresponding hyperon fractions are smaller. 



III. THE BULK VISCOSITY 



In general bulk viscosity is a dissipative process in which the compression is converted into 
heat. It is due to the instantaneous difference between the local physical pressure p and the 
thermodynamic pressure p 

p — p = — C V.-y, 

where v is the local velocity. In neutron star matter one is interested in describing small deviations 
(oscillations) 6xi = Xi — Xi around the thermodynamic equilibrium, characterized by the variables 
(densities)ij. 

We will now turn to the calculation of C, for the specific case of a neutron star with hyperons. 
Since the three reactions involve neutrons the fiuctuation of the neutron fraction x„ can be 
used as the primary parameter. For this situation a general expression for the bulk viscosity has 
been derived in Ref. ^ 

C = :-^(|^) ^, (12) 

with T the relaxation time, uj the pulsation frequency of the neutron star, p the pressure, ns the 
baryon density. 

We restrict ourselves to the densities, at which nucleons, S~'s, and A's are present. Defining 
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the chemical imbalance as 



X = 2 l^n - f^p - Ht.- = f^n - fJ'A, (13) 

and assuming that in first order the difference AF^ between the rates for the various direct reactions 
Ta and inverse reactions Ta in Eqs. is proportional to x 

AT = J2 ATa = X X, (14) 

a 

with A determining the viscosity. Hence, the relaxation time r can be expressed in terms of the 
variation of the imbalance with the neutron density fluctuation 

UBX ^^^^ 



AT (dx/dxn) \{dx/dxn)' 
Thus the main task is to evaluate the dx/dxn- We wish to include several weak nonleptonic hyperon 
reactions. To do this in a proper way one has to take into account that the four relevant baryon 
number densities, XmXp,x\^ and xs, are related to each other by the following three constraints: 
baryon conservation, charge neutrality, and chemical equilibrium for strong processes and in par- 
ticular for the reaction n + A<-^p + S^. By using these conditions one can express x iii terms of 
the fractions Xj and assuming that all leptonic reaction rates are much smaller than those which 
produce the hyperon bulk viscosity, one obtains in the density region with npT,^ matter 

dx 

= 2 ann - ( C^pn + Q^E-n + <^np + ^nE" ) 

CLXfi 

+ (1/2) {app + "s-p + "pS- + "s-s-), (16) 
and in the region with npT,~A matter 

dx ^ ^ _^ { Pn - Pa )( anp - QAp + QnS- " "AS-) 

dxn 2 (3a - Pp - /3s- 

{2 Pn - Pp - Pt.- )( "nA - "AA ) ..„v 

- 2P, - P, - p^. ' 

where the a's correspond to the partial derivatives of the chemical potentials with respect to the 
various fractions 

"ii = (18) 

and 

Pi = oini + oiM - api - as-i) (19) 
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where i and j stand for n, p, 'E , and A. These quantities are obtained from the EoS. 
The real part of the hyperon bulk viscosity 

1 + 7-2 V dx.^ J duB^ 



(20) 



where r is given in Eq. (|15j) . The core of the neutron star is assumed to pulsate with a frequency of 
about uj ~ 10^ — 10^ s^^. In the high frequency limit 1 << ujt the bulk viscosity Q is proportional 
to the inverse of the relaxation time, C, oc r^^ oc A, whereas in the low frequency limit 1 >> wr 
the bulk viscosity is proportional to the relaxation time, oc r oc A~^. In the following section the 
difference in the rates given by AT in Eq. ()14() in first order in x is derived, from which A = AF/x 
can be obtained. 

IV. COLLISION RATE 

In this section the various rates AFjyi^jjv'y' foi^ the processes A^i + A''2 N' + Y' are considered, 
where the is a nucleon and the y is a hyperon. First we consider the one-pion exchange in Born 
approximation. For completeness the rates using the contact W exchange interaction jQ] are given 
in section IIV Bl 

A. One-pion exchange 
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FIG. 3: The pion exchange diagrams in Born approximation for the weak nonleptonic hyperon processes 
with only one hyperon involved. 



In order to obtain AFjvi]V2jv'y'> first the rate ^ nin2-.n'y' is considered. The relevant free space 
diagrams are shown in Fig. |31 The strangeness changing weak vertex is given by 

F^^, = Gpmli A^y. + B^y, 75 ), (21) 

for which the constants A and B determine the strengths of the parity violating and parity con- 
serving y — > A^ + TT amplitudes, respectively, and are specified in section The strong vertex 
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FIG. 4: The two types of closed diagrams at the two-loop level. 



-^j^jv' — S'jviv' 75 



(22) 



To compute the rates in the medium one needs to account for Pauli blocking; it is convenient 
to use the optical theorem to convert the free space diagrams of Fig. El to the closed diagrams of 
Fig. ^ Using finite temperature Green functions in the QPA, one can express the collision rate as 

1 f d^p^, d'^p^, d^p^, dW 



]VlJV2-.JV'i" 



S J (2 7r)4 (2^)4 (2 7r)4 {2 
Z {2ttY 6^{pni +Pn2 -Py' -Pn')- 



(23) 



In Eq. (ESI) one has 



(24) 



with 



TV 



Gjvi^(Pivi) -^^j^ijv' G~^> {Pn') Fnin' 



^ N2Y' ^y' \Py' ) ^ N2Y' 



(25) 



with DT^[k'^) = ^ , and ki = pjv, — pj^i (for i = 1, 2), corresponding to the diagram on the left 



in Fig. 0]and 



Tr 



D^{kl) D^{kl) + {A^i ^ N2} 



(26) 



corresponding to that of the right in Fig. The symmetry factor S is for the pnpA process 5 = 1 
and for the nnpTi and nnnA processes S = 2. After evaluation of the p^ integrals in the QPA, one 
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obtains 

(fpNi d?PN2 cP'Pn' d^Py' 



N1N2^n'y' — J 



2 mjvi (2 7r)3 2 m^va (2 n)^ 2 m^v' (2 7r)3 2 my/ (2 7r)3 



^ \MniN2N'y\ Ini /]V2 (1 /jv') (1 Iy') 
(2 + - - 

{Pni + Pn2 - Pn' - Py' ) , (27) 

where pi and Ei are the particle momentum and energy, respectively, and fi = {1 + exp[(ei — 
fj.i)/T]}^^ is the Fermi-Dirac distribution function. The matrix element is 

-^N1N2N'y' ~ [ ''^JV' -^JVliv' ^JVl "Uy' Fjsi2y' '^N2 -^^77(^1) 
^N2N' "^N2 Uyl F]siiy' 

Uni ]. (28) 

The evaluation of rjvi]V2^]v'y' with the standard technique takes advantage of the strong de- 
generacy of the participating particles in neutron star matter. The multidimensional integral is 
decomposed into an energy and an angular integration. All momenta pi are placed on the appro- 
priate Fermi spheres wherever possible. Furthermore, we introduce the dimensionless quantities 

= ^^^;C = f (29) 

The rate in Ea. ()27() can be factorized in the form 



with 

m 

and 



^N1N2^n'y' — '^[Jin2N'y' -^(0) ("^0) 



^ r'=° 4 TT^ f + f3 

n / dmfiy^)SiJ2y^ + = (31) 

i=i -^-"^ i=i e'- - I D 



p(o) _ PFN. T3 f 9(1 - \C{9^>,e,,)\) 

^n^n2n'y' - 8{2 7rfJ ^'^-'^'^^^ ^1 _ CH0N',0Y') 

^ X] \^N1N2N'y'\'^, (32) 
spin 

where pfni is the Fermi momentum of baryon Ni, O^i is the angle between p^ri and p^i ^ By' is the 
angle between p^vi and py' , and 



C{e^>,9y>) = {-Pfn^ - Pfj^> - Ppy' + 2 pfni Pfn' cos( 



'JV' 



+ 2 pfni Pfy' cos{9y') - 2 pp^f ppy, cos{9 JV'j COS(f7y/j 

+ pI^,)/{2pI^, sin(^^0 sin(0yO)- (33) 
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The conditions in section^ Eqs. apply in thermodynamic equihbrium. However, small 

deviations from thermodynamic equilibrium occur just after a neutron star is born. For the 
bulk viscosity the small deviations from the weak chemical equilibrium are important. For this 
purpose we need to consider AF^i^jiv'y' = rjvijv2_jv'y' — r]v'i'-.ivijv2i where the inverse rate is 
'^n'y^nin2 = ^^nin2n'y' Thcrcfore, in the linear approximation, valid for small deviations 

1^1 = \x\/T « 1 one has 

with 

AI = m - H-O = (35) 

The AI is the same for the weak nonleptonic hyperon processes of Eqs. H3I5() . The total rate 
considered for the hyperon bulk viscosity is 

AT = AFpnpA + AFnnnA + 2 Ar^npT,-- (36) 

This result is inserted into Eq. H15|) to obtain the relaxation time. 



B. Contact W exchange interaction 




FIG. 5: The diagram of the contact W exchange interaction for the nnpY. process and the pnpA process . 

In order to illustrate the difference between the pion exchange approach and the contact W 
exchange picture, we also consider the latter. The rates for W exchange have been calculated 
before by Haensel et al. P| and by Lindblom and Owen however in the angle-averaged approx- 
imation. Performing an analogous derivation as in the previous section, one obtains the expression 
in Eq. (1151) . For convenience of the reader here we give the simpler matrix elements in the nonrel- 

n 

ativistic limit (for the general expression for the matrix elements we refer to Ref. |6() 
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WnnpJi-? = 8 Gl sm2(2 Oc) ml m^- ( 1 + 3 cj cf'~f (37) 

and 

\MpnpA.? = 8 Gl sin2(2 0c) mA ( 1 + 3 jcjp ), (38) 

respectively. Note that no meson propagators D{k'^) are present in Eqs. (|37p and (|38() . i.e. it has 
a different density dependence in the medium than the OPE matrix elements. The nnnA process 
has no simple W exchange contribution. The results are shown in the following section. 



V. RESULTS AND DISCUSSION 



Some remarks about the differences between the meson exchange and the W exchange picture 
are appropriate. First in the OPE (in which phenomenological input is used for both the weak and 
strong ttNN couplings) the pion exchange leads to a finite range interaction and the pseudo scalar 
ttNN coupling to the presence of momentum dependent terms, pi/m. These effects are absent in 
the simple W point coupling. Hence, in addition to a different overall strength one also expects to 
find a different density dependence of the matrix elements for the two processes in the medium. 
Secondly, the neutral current process, such as nnnA process, cannot be described in the simple W 
exchange approach. 

The EoS applied for the following figures is given by Balberg and Gal [3| with parametrization 
7 = (5 = 2 (case 1) and also j = 6 = 5/3 (case 2). In the calculation for the case of pion exchange 
we use the following values for the strong and weak vertices: 5^„/(47r) = 14.2 and gnp = V^gnn- 
The values A^^- = 1.93, A^a = -1.07, ApA = 1.46, B^^- = -0.63, 4iA = -7.19, and BpA = 9.95 
are based upon the experimental data on pionic decay of hyperons [isj l [l^ . In the W exchange 
picture the standard values of the weak coupling constants are = —qa = —1.26, c^^ = —0.72, 
and = 0.34 and Gf is the Fermi weak coupling constant. These values are also used by 

• Q n 

Lindblom and Owen [g], while Haensel et a/Jj] treats the interaction matrix as a free parameter. 



Compared the work of Lindblom and Owen [g] and Haensel et al. |J] effects of superfluidity are not 
taken into account and there is no angle averaging in this work. 

In Fig. El we compare the rates for the n + n ^ p + S~ process for the W exchange and the 
pion exchange as a function of density for case 1 of the EoS. Concerning the W exchange, we 
note that the nonrelativistic result is a poor approximation to the relativistic result since there is 
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FIG. 6: The equilibrium rate T for n + n^p + S as a function of the baryon number density at T = 10 K. 



an almost complete cancellation between the leading order vector and axial- vector contributions, 
(1 + Sc^c^ )^ ~ 1/16. The overall rate calculated with OPE is almost two orders of magnitude 
larger than the one corresponding to relativistic W exchange. The difference can be attributed to 
larger values of the couplings, Any' , Bny' , and qnn' > 1, in the OPE case. At low momentum 
transfer (low density), this effect would be compensated for by p-wave character of the ttNN 
coupling, a.k. However, for momenta relevant in neutron stars this does not happen. In addition, 
a stronger density dependence is observed, which is related to finite range pion exchange in Eq. (|28() . 

Because of the rapid increase of the A fraction after its appearance, the proton and fraction 
drop abruptly and also their Fermi momenta. Therefore, the kink in the rate can be attributed to 
the appearance of the A at ub ~ 0.33 fm~^. 

The rates for pnpA process are shown in Fig. [Tj As mentioned above the OPE shows a stronger 
density dependence because of the presence of the meson propagator. Although at baryon densities 
where the hyperons are present a partial cancellation between the direct and exchange OPE matrix 
elements occurs, the overall OPE rate is about one order of magnitude larger than the one obtained 
with W exchange, which again can be attributed to the larger values of the meson couplings. 

As mentioned above the neutral current nnnA process does not receive a contribution from W 
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FIG. 7: The equilibrium rate T for p + n ^ p + A and n + n^n + Aasa function of the baryon number 
density at T = 10^ K. 



exchange. In the OPE approach the magnitude of this rate is a factor 2-3 smaher than that of the 
pnpA process. We note that there exists experimental information on the ratio of the nnnA and 
pnpA processes in the weak decay of hypernuclei which suggests that the simple OPE mechanism 



The inclusion of other mesons, such as kaons, could improve the 



is not sufficiently accurate 
situation. 

The relaxation time is shown as a function of the baryon density in Fig. |H1 One sees that the 
relaxation time in the OPE picture is about two orders of magnitude smaller than in the W picture. 
In the latter case after the appearance of S~, the relaxation time is determined by the nnpT,~ 
process; after the occurrence of the A, the pnpA process takes over and it dominates the relaxation 
time. Therefore, a drop in the relaxation time occurs at the appearance density of the A. In the 
OPE picture after the appearance of the A, the nnpT,^ process remains the most important process 
for the relaxation time. Therefore at the appearance density of the A a small abrupt increase in 
the relaxation time occurs. 

To compute the viscosity one needs to assume a value for the frequency; typical values of the 
frequency of the pulsations are uj = 10^ - 10^ s-^. We have used to = 10^ s'^ in Fig. H For the 
OPE picture, the high frequency limit is applicable. Whereas for W exchange at T = 10^ K and 
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FIG. 8: The relaxation time r as a function of the baryon number density at T = 10 K. 
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FIG. 9: The bulk viscosity C as a function of the baryon number density. 
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FIG. 10: The bulk viscosity C as a function of the baryon number density for different equations of state at 
T = IQi" K. 



at low density we are even in the low frequency limit. The bulk viscosity in the OPE picture is 
about 1-2 orders of magnitude smaller than that in the W exchange picture. 

To investigate the sensitivity of the bulk viscosity to details of the EoS we also considered 
different values of the parameters. As can be seen in Fig. EH the bulk viscosity is rather insensitive 
to the details of the EoS except with respect to the appearance of the A. 



VI. CONCLUSION 



In this paper the bulk viscosity due to weak nonleptonic hyperon processes has been studied. 
This viscosity is relevant in connection with damping of neutron star pulsations, especially in the 
damping of r-modes. In particular we considered pion exchange, which is based upon empirical 
input, in Born approximation to describe the weak nonleptonic hyperon processes instead of the 
W exchange Q, Q, 0]. The conclusions are: i) The bulk viscosity in the OPE picture is about 1-2 
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orders of magnitude smaller than that in the W exchange picture, ii) The rates of the nnp'E and 
pnpA process are of the same order of magnitude using OPE. This result is in contrast with the 
case of W exchange, because the pnpA process is one order of magnitude larger than nnpT, process, 
iii) The nnnA process can be included in the calculation for the bulk viscosity using OPE. iv) The 
bulk viscosity is rather insensitive for the EoS used except with respect to the appearance of the 
A. 

Therefore, the contact W exchange interaction is probably to naive for quantitative calculations. 
In a more realistic OPE approach the contributions of the hyperons to the bulk viscosity are less 
pronounced. 
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